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Abstract 

Let S' be a planar n-point set. A triangulation for S is a maximal plane straight-line graph 
with vertex set S. The Voronoi diagram for S is the subdivision of the plane into cells such 
that all points in a cell have the same nearest neighbor in S. Classically, both structures can 
be computed in O(nlogn) time and 0{n) space. We study the situation when the available 
workspace is limited: given a parameter s S n}, an s-workspace algorithm has read¬ 

only access to an input array with the points from S in arbitrary order, and it may use only 
0(s) additional words of 0(logn) bits for reading and writing intermediate data. The output 
should then be written to a write-only structure. We describe a deterministic s-workspace 
algorithm for computing an arbitrary triangulation of S in time 0{n^js + nlognlogs) and 
a randomized s-workspace algorithm for finding the Voronoi diagram of S in expected time 
0((n^/s) logs -I- nlogslog* s). 


1 Introduction 

Since the early days of computer science, a major concern has been to cope with strong memory 
constraints. This started in the ’70s |22j when memory was expensive. Nowadays, a major motiva¬ 
tion comes from a proliferation of small embedded devices where large memory is neither feasible 
nor desirable (e.g., due to constraints on budget, power, size, or simply to make the device less 
attractive to thieves). 

Even when memory size is not an issue, we might want to limit the number of write opera¬ 
tions: one can read flash memory quickly, but writing (or even reordering) data is slow and may 
reduce the lifetime of the storage system; write-access to removable memory may be limited for 
technical or security reasons (e.g., when using read-only media such as DVDs or to prevent leaking 
information about the algorithm). Similar problems occur when concurrent algorithms access data 
simultaneously. A natural way to address this is to consider algorithms that do not modify the 
input. 

*WS and PS were supported in part by DFG Grants MU 3501/1 and MU 3501/2. YS was supported by the DFG 
within the research training group “Methods for Discrete Structures” (GRK 1408). 
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The exact setting may vary, but there is a common theme: the input resides in read-only memory, 
the output must be written to a write-only structure, and we can use 0(s) additional variables to 
find the solution (for a parameter s). The goal is to design algorithms whose running time decreases 
as s grows, giving a time-space trade-off |23) . One of the first problems considered in this model 
is sorting mm- Here, the time-space product is known to be H(n^) [8], and matching upper 
bounds for the case b G H(logn) H 0(n/ log n) were obtained by Pagter and Rauhe |21| (6 denotes 
the available workspace in bits). 

Our current notion of memory constrained algorithms was introduced to computational ge¬ 
ometry by Asano et al. j^, who showed how to compute many classic geometric structures with 
0(1) workspace (related models were studied before |9]). Later, time-space trade-offs were given 
for problems on simple polygons, e.g., shortest paths |1], visibility |6], or the convex hull of the 
vertices [5]. 

We consider a model in which the set S' of n points is in an array such that random access to 
each input point is possible, but we may not change or even reorder the input. Additionally, we have 
0{s) variables (for a parameter s G {1,..., n}). We assume that each variable or pointer contains a 
data word of 0(logn) bits. Other than this, the model allows the usual word RAM operations. In 
this setting we study two problems: computing an arbitrary triangulation for S and computing the 
Voronoi diagram VD(S) for S. Since the output cannot be stored explicitly, the goal is to report 
the edges of the triangulation or the vertices of VD(S) successively, in no particular order. Dually, 
the latter goal may be phrased in terms of Delaunay triangulations. We focus on Voronoi diagrams, 
as they lead to a more natural presentation. 

Both problems can be solved in O(n^) time with 0(1) workspace |1] or in 0(n log re) time with 
0(re) workspace [7|- However, to the best of our knowledge, no trade-offs were known before. Our 
triangulation algorithm achieves a running time of 0{n^ / s -\- re log re logs) using 0(s) variables. 
A key ingredient is the recent time-space trade-off by Asano and Kirkpatrick for triangulating a 
special type of simple polygons jSj. This also lets us obtain significantly better running times for 
the case that the input is sorted in x-order; see Section For Voronoi diagrams, we use random 
sampling to find the result in expected time 0((re^logs)/s-|-relogslog* s)); see Section]^ Together 
with recent work of Har-Peled m, this appears to be one of the first uses of random sampling 
to obtain space-time trade-offs for geometric algorithms. The sorting lower bounds also apply to 
triangulations and Voronoi diagrams (since we can reduce the former to the latter). This implies 
that our second algorithm is almost optimal. 

2 A Time-Space Trade-Off for General Triangulations 

In this section we describe an algorithm that outputs the edges of a triangulation for a given point 
set S in arbitrary order. For ease in the presentation we first assume that S is presented in sorted 
order. In this case, a time-space trade-off follows quite readily from known results. We then show 
how to generalize this for arbitrary inputs, which requires a careful adaptation of the existing data 
structures. 

2.1 Sorted Input 

Suppose the input points S = {g^i,..., Qn} are stored in increasing order of x-coordinate and that 
all x-coordinates are distinct, i.e., Xi < Xj+i for I < f < re, where Xi denotes the x-coordinate of qi. 
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A crucial ingredient in our algorithm is a recent result by Asano and Kirkpatrick for triangulating 
monotone mountain^ (or mountains for short). A mountain is a simple polygon with vertex 
sequence vi,V 2 , ■ ■ ■ ,Vk such that the x-coordinates of the vertices increase monotonically. The edge 
viVk is called the base. Mountains can be triangulated very efficiently with bounded workspace. 

Theorem 2.1 (Lemma 3 in |3], rephrased). Let H be a mountain with n vertices, stored in sorted 
x-order in read-only memory. Let s G {2, ... ,n}. We can report the edges of a triangulation of H 
in O(nlog^n) time and using 0{s) words of space. 

Since S is given in x-order, the edges for 1 < i < n, form a monotone simple polygonal 

chain. Let Part (S') be the subdivision obtained by the union of this chain with the edges of the 
convex hull of S (denoted by conv(S)). A convex hull edge is long if the difference between its 
indices is at least two (i.e., the endpoints are not consecutive). The following lemma (illustrated in 
Fig. 0 lets us decompose the problem into smaller pieces. 



Figure 1: Any face of Part(S) is a mountain that is uniquely associated with a long convex hull 
edge. 


Lemma 2.2. Any bounded face of Part(S) is a mountain whose base is a long convex hull edge. 
Moreover, no point of S lies in more than four faces o/Part(S). 

Proof. Any point Qi € S has at most four neighbors in Part(S): Qi-i, Qi+i, its predecessor and its 
successor along the convex hull (if qi lies on conv(5)). Thus, no point of S belongs to more than 
four faces of Part (5). 

Next we show that every face F of Part(S') is a mountain with a long convex-hull edge as its 
base. The boundary of F contains at least one long convex-hull edge e = {qi,qj) {i < j), as other 
edges connect only consecutive vertices. Since the monotone path qi,... ,qj forms a cycle with the 
edge e and since the boundary of T is a simple polygon, we conclude that e is the only long convex- 
hull edge bounding F. Recall that e is a convex hull edge, and thus all points qi+i, ..., qj-i lie on 
one side of e and form a monotone chain (and in particular T is a mountain with base e). □ 

The algorithm for sorted input is now very simple. We compute the edges of the convex hull 
(starting from the leftmost point and proceeding in clockwise order). Whenever a long edge would 
be reported, we pause the convex hull algorithm, and we triangulate the corresponding mountain. 
Once the mountain has been triangulated, we resume with the convex hull algorithm until all convex 

^Also known as unimonotone polygons |15| . 
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hull edges have been computed. The trade-off now follows from already existing trade-offs in the 
various subroutines. 


Theorem 2.3. Let S be a set of n points, sorted in x-order. We ean report the edges of a triangu¬ 
lation of S in 0{n^) time using 0(1) variables, in 0(n^ log n/2'^) time using 0{s) variables (for any 
s G 0(loglogn)no(logn)J; and in O(nlogpn) time using 0{p\ogpn) variables (for any 2 < p < n). 


Proof. Correctness follows from Lemma 2.2, so we focus on the performance analysis. The main 
steps are: (i) computing the convex hull of a point set given in x-order; and (ii) triangulating a 
mountain. 

By Theorem 2.1 we can triangulate a mountain Fi with m vertices in time 0(ni log^ n*) with 
0{s) variables. We do not need to store T) explicitly, since its vertices constitute a consecutive 
subsequence of S and can be specified by the two endpoints of the base. No vertex appears in 
more than four mountains by Lemma |2.2[ so the total time for triangulating the mountains is 

0{ni log^ Ui) = 0{nlogg n). By reusing space, we can ensure that the total space requirement is 
0{s). 

Now we bound the time for computing conv(S'). This algorithm is paused to triangulate moun¬ 
tains, but overall it is executed only once. There are several convex hull algorithms for sorted point 
sets under memory constraints. If s G ©(1), we can use gift-wrapping (Jarvis march |l7)L which 
runs in 0(n?) time. Barba et al. |5] provided a different algorithm that runs in 0(n^ logn/2'^) time 
using 0{s) variables (for any s G o(logn))j^ This approach is desirable for s G ll(loglogn)no(logn). 
As soon as s = n(logn), we can use the approach of Chan and Chen |in| . This algorithm runs in 
0(n logpu) time and uses O(plogpn) variables, for any 2 < p < n. Regardless of the size of the 
workspace, the time for computing the convex hull dominates the time needed for triangulating all 
mountains. □ 


A similar approach is unlikely to work for the Delaunay triangulation, since knowing the x-order 
of the input does not help in computing it mi. 


2.2 General Input 


The algorithm from Section 2.1 uses the sorted order in two ways. Firstly, the convex-hull algorithms 
of Barba et al. |5] and of Chan and Chen |10) work only for simple polygons (e.g., for sorted input). 
Instead, we use the algorithm by Darwish and Elmasry [13] that gives the upper (or lower) convex 
hull of any sequence of n points in 0(n^/(s logn) -|- nlogn) time with 0(s) variable^ matching 
known lower bounds. Secondly, and more importantly, the Asano-Kirkpatrick (AK) algorithm for 
triangulating a mountain requires the input to be sorted. To address this issue, we simulate sorted 
input using multiple heap structures. This requires a close examination of how the AK-algorithm 
accesses its input. 

Let T be a mountain with n vertices. Let and F-^ denote the vertices of F in ascending and 
in descending x-order. The AK-algorithm has two phases, one focused on and the other one on 


^In fact, Barba et al. show how to compute the convex hull of a simple polygon, but also show that both problems 
are equivalent. The monotone chain can be completed to a polygon by adding a vertex with a very high or low 
j/-coordinate. 

^Darwish and Elmasry 1131 state a running time of 0(n^/s + nlogn), but they measure workspace in bits, while 
we use words. 
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Each pass computes a portion of the triangulation edges, uses 0{s) variables, and scans the 
input 0(log^n) times. We focus on the approach for F^. 

As mentioned, the algorithm uses 0(loggn) rounds. In round i, it partitions F into blocks of 
0{\F\/s^) consecutive points that are processed from left to right. Each block is further subdivided 
into 0{s) sub-blocks bi,...,bk of size The algorithm does two scans over the sub¬ 

blocks. The first scan processes the elements from left to right. Whenever the first scan finishes 
reading a sub-block 5j, the algorithm makes active and creates a pointer li to the rightmost 
element of bt. The second scan goes from right to left and is concurrent to the first scan. In each 
step, it reads the element at li in the rightmost active sub-block 6j, and it decreases U by one. If li 
leaves bi, then bi becomes inactive. As the first scan creates new active sub-blocks as it proceeds, 
the second scan may jump between sub-blocks. The interested reader may find a more detailed 
description in 

To provide the input for the AK-algorithm, we need the heap by Asano et al. [2]. For complete¬ 
ness, we briefly restate its properties here. 

Lemma 2.4 (|2]). Let S be a set of n points. There is a heap that supports insert and extract-min 
(resp. extract-max) in 0((n/(slogn) -|-logs)I?(n)) time using 0{s) variables, where D{n) is the 
time to decide whether a given element currently resides in the heap (is alivedj^ 

Proof. We first describe the data structure. Then we discuss how to perform insertions and extract- 
min operations. 

We partition the input into slogn consecutive buckets of equal size, and we build a complete 
binary tree T over the buckets. Let u be a node of T with height h. Then, there are 2^ buckets 
below V in T. We store 2h information bits in v to specify the minimum alive element below v. 
The first h bits identify the bucket containing the minimum. We further divide this bucket into 
2^ consecutive parts of equal size, called quantiles. The second h bits in v specify the quantile 
containing the minimum. If 2h > logn, we use logn bits to specify the minimum directly. Hence, 
the total number of bits is bounded by 

min{2/i, log n} = O(slogn). 

h=0 


log(slogn) 


E 


s log 
2 >^ 


Therefore we need 0{s) variables in total. 

Let u be a node with height h. To find the minimum alive element in T below v, we use the 2h 
information bits stored in v. First, we identify the bucket containing the minimum and the correct 
quantile within this bucket. This quantile contains elements. For each element in the 

quantile, we decide in D(n) time whether it is alive, and we return the minimum such element. This 
takes 0( ^fa^”^g^ Zj(n)) time in total. 


insert: Assume we want to insert an element x that is at position i in the input array. Let v be the 
parent of the leaf of T corresponding to the bucket that contains x. We update the information 
bits at each node u on the root path starting at u. To do so, we use the information bits in u to 


reduce triangulation to the next smaller right neighbor (NSR) and the next smaller left neighbor (NSL) 
problem and present an algorithm for NSR if the input is in a;-order. This implies an NSL-algorithm by reading the 
input in reverse. 

®The bounds in [5] do not include the factor D{n) since the authors studied a setting similar to Lemma 
it takes 0(1) time to decide whether an element is alive. 


2.5 


where 
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find the minimum element in the buckets covered by u, as described above. Then we compare 
it with X. If X is larger, we are done and we stop the insertion. Otherwise, we update the 
information bits at u to the bucket and quantile that contain x. If we reach and update the 
root node, we also update the pointer that points to the minimum element in the heap. The 
work per node is dominated by the costs for finding the minimum, which is 
Thus, the total cost for insertion is bounded by 

log(slogn) 

Y -ITT - D{n)=o(Y^D{n)). 

^ 2^slogn Vslogn / 

h=o ^ ^ 

extract-min: First we use the pointer to the minimum alive element to determine the element x 
to return. Then we use a similar update strategy as for insertions. Let v be the leaf node 
corresponding to the bucket of x. We first update the information bits of v by scanning through 
the whole bucket of v and determining the smallest alive element. Since a bucket contains 
0(n/s log n) elements, this needs time O{n/{s log n)D(n)). Then we update the information 
bits of each node u on the path for v as follows: let vi and t >2 be the two children of u. We 
determine the minimum alive element in the buckets covered by and V 2 , take the smaller 
one, and use it to update the information bits at u. Once we reach the root, we also update 
the pointer to the minimum element of the heap to the new minimum element of the root. 
The total time again is bounded by 


□ 

Lemma 2.5 (|2]). Let S be a set of n points. We can build a heap with all elements in S in 0{n) 
time that supports extract-min in 0{n/{slogn) +logn) time using 0(s) variables. 

Proof. The construction time is given in [2]. To decide in 0(1) time if some x € S is alive, we store 
the last extracted minimum m and test whether x > m. □ 


We now present the complete algorithm. We show how to subdivide S into mountains T) and 
how to run the AK-algorithm on each FI. By reversing the order, the same discussion applies to Ff. 
Sorted input is emulated by two heaps Lfi, H 2 for S according to x-order. By Lemma [2. 5 [ each heap 
uses 0{s) space, can be constructed in 0(n) time, and supports extract-min in 0(n/(s log n) -|-log re) 
worst-case time. We will use Hi to determine the size of the next mountain T) and H 2 to process 
the points of Fi. 

We execute the convex hull algorithm with 0(s) space until it reports the next convex hull edge 
pq. Throughout the execution of the algorithm, heaps Hi and H 2 contain exactly the points to the 
right of p. We repeatedly extract the minimum of Hi until q becomes the minimum element. Let 
k be the number of removed points. 

If /c = 1, then pq is short. We extract the minimum of H 2 , and we continue with the convex 
hull algorithm. If /c > 2, then Lemma 2.2 shows that pq is the base of a mountain F that consists 
of all points between p and q. These are exactly the k -\- 1 smallest elements in H 2 (including p 
and q). If /c < s, we extract them from H 2 , and we triangulate F in memory. If A; > s, we execute 
the AK-algorithm on F using 0{s) variables. At the beginning of the ith round, we create a copy 
of H 2 , i.e., we duplicate the 0{s) variables that determine the state of H 2 . Further, we create 
an empty max-heap iL(ij) using 0(s) variables to provide input for the second scan. To be able 
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to reread a sub-block, we create a further copy of H 2 . Whenever the AK-algorithm requests 
the next point in the first scan, we simply extract the minimum of When a sub-block is fully 
read, we use to reread the elements and insert them into -ff(jj). Now, the rightmost element of 
all active sub-blocks corresponds exactly to the maximum of One step in the second scan is 

equivalent to an extract-max on Hyy. 

At the end of a round, we delete Hyy and Hyy, so that the space can be reused in the 
next round. Once the AK-algorithm finishes, we repeatedly extract the minimum of H 2 until we 
reach q. 


Theorem 2.6. We can report the edges of a triangulation of a set S of n points in time 0{n^/s + 
n log re logs) using 0{s) additional variables. 


Proof. Similarly as before, correctness directly follows from Lemma 2.2 and the correctness of the 
AK-algorithm. The bound on the space usage is immediate. 

Computing the convex hull now needs 0(re^/(s log re) -|- re log re) time |13) . By Lemma 2.5 the 
heaps Hi and H 2 can be constructed in 0(re) time. During execution, we perform re extract-min 
operations on each heap, requiring 0(re^/(s log re) -|- re log re) time in total. 


Let 


Fj be a mountain with Uj vertices that is discovered by the convex hull algorithm. If 


re 


j < s, then Fj is triangulated in memory in 0{nj) time, and the total time for such mountains 
is 0{n). If Uj > s, then the AK-algorithm runs in 0(nj log^ nj) time. We must also account 
for providing the input for the algorithm. For this, consider some round i > 1. We copy H 2 to 
in 0{s) time. This time can be charged to the first scan, since nj > s. Furthermore, we 
perform Uj extract-min operations on Hyy Hence the total time to provide input for the first scan 
is 0(rejre/(s logre)-|-rej log re). 

For the second scan, we create another copy of H 2 . Again, the time for this can be charged to 
the scan. Also, we perform nj extract-min operations on which takes 0{njnj(s log re) -|-rej log re) 
time. Additionally, we insert each fully-read block into Hyy. The main problem is to determine if an 
element in Hyy is alive: there are at most 0{s) active sub-blocks. For each active sub-block 5,, we 
know the first element yi and the element Zi that Zj points to. An element is alive if and only if it is 


in the interval [?/j, Zi\ for some active This can be checked in 0(log s) time. Thus, by Lemma 2.4 
each insert and extract-max on Hyy takes 0((re/(slogre) -|-logs) logs) time. Since each element is 
inserted once, the total time to provide input to the second scan is 0(rej log(s)(re/(s logre) -|-logs)). 
This dominates the time for the first scan. There are ©(log^rej) rounds, so we can triangulate 
Fj in time OiynjloggUj + nj\og{nj){n/{slogn) -|-logs)). Summing over all Fj, the total time is 


0{n^ j s -|- re log re log s) 


□ 


3 Voronoi Diagrams 

Given a planar re-point set S', we would like to find the vertices of VD(S). Let K = {pi,P2,P3} be a 
triangle with SniL = 0, S C conv(A'), and so that all vertices of VD(S) are vertices of VD(SUiL). 
For example, we can set K = {(—re, — k), (— k, re), (0, re)} for some large re > 0. Since the desired 
properties hold for all large enough re, we do not need to find an explicit value for it. Instead, 
whenever we want to evaluate a predicate involving points from K, we can take the result obtained 
for re —>• 00 . 

Our algorithm relies on random sampling. First, we show how to take a random sample from 
S with small workspace. One of many possible approaches is the following one that ensures a 
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worst-case guarantee: 


Lemma 3.1. We can sample a uniform random subset R <Z S of size s in time 0{n + slogs) and 
space 0{s). 

Proof. The sampling algorithm consists of two phases. In the hrst phase, we sample a random 
sequence / of s distinct numbers from [n]j^ The phase proceeds in s rounds. At the beginning of 
round k, for k = 1,..., s, we have already sampled a sequence I of A; — 1 numbers from [n], and we 
would like to pick an element from [n] \ I uniformly at random. We store I in a binary search tree 
T. We maintain the invariant that T stores with each element xG[n — A:-|-l]n/a replacement 
Px £ {n — A; -|- 2,..., n} \ I such that [n] \ / = ([n — A: -|- 1] \ /) U {px | x G [n — A; -|- 1] H /}, see 
Figure 1^ In round k, we sample a random number x from [n — A; -1- 1], and we check in T whether 



Figure 2: Sampling a random sequence I from [n]. At the beginning of round k, we have already 
sampled A: — 1 elements (shown in gray). Each element x G [n — A:-|-l]n / has a replacement 
Px G {n — A: -|- 2,..., n} \ I (indicated by the arrows). In round k, we pick a random number 
X G [n — A: -|- 1]. If X is already contained in I, we add px to I. Otherwise, we add x. 


X G /. If not, we add x to I (and T), otherwise, we add px to I (and T). By the invariant, we add 
a uniform random element from [re] \ / to /. 

It remains to update the replacements, see Figure If x = re — A: -|- 1, we do not need a 
replacement for x. Now suppose x<re — A:-|-l. Ifre — A:-|-10/, we set pa, = re — A: -|- 1. Otherwise, 
we set Px = Pn-k+i- This ensures that the invariant holds at the beginning of round A: -|- 1, and it 
takes O(logs) time and 0{s) space. We continue for s rounds. At the end of the hrst phase, any 
sequence of s distinct numbers in [re] is sampled with equal probability. Furthermore, the phase 
takes 0(slogs) time and 0(s) space. 

In the second phase, we scan through S to obtain the elements whose positions correspond to 
the numbers in I. This requires 0{n) time and 0{s) space. □ 


We use Lemma 3.1 to hnd a random sample R Q S of size s. We compute VD(i?UiF), triangulate 
the bounded cells and construct a planar point location structure for the triangulation. This takes 
O(slogs) time and 0(s) space jl^. By our choice of K, all Voronoi cells for points in R are 
bounded, and every point in S lies in a bounded Voronoi cell. Given a vertex v G VD(i? U K), the 
conflict circle of v is the largest circle with center v and no point from ii U AT in its interior. The 
conflict set of v contains all points from S that he in the conhict circle of v, and the conflict 
size by of V is \By\. We scan through S to hnd the conhict size by for each vertex v G VD(i? U K): 
every Voronoi vertex has a counter that is initially 0. For each p G S \ (RU K), we use the point 
location structure to hnd the triangle A of VD(AU K) that contains it. At least one vertex v of 
A is in conhict with p. Starting from v, we walk along the edges of VD(i? U K) to hnd all Voronoi 
vertices in conhict with p (recall that these vertices induce a connected component in VD(i?U K)). 
We increment the counters of all these vertices. This may take a long time in the worst case, so we 


We write [n] for the set {1,..., n}. 
































Figure 3: Finding a replacement for x. If x = n — fc + 1, we do not need a replacement for x in the 
next round (top left). If n — /c + 1 is not sampled yet, we can make it the replacement for x (top 
right). Otherwise, we make the old replacement for n — k + 1 the new replacement for x (bottom). 


impose an upper bound on the total work. For this, we choose a threshold M. When the sum of the 
conflict counters exceeds M, we start over with a new sample R. The total time for one attempt 
is 0{nlogs + M), and below we prove that for M = 0(n), the success probability is at least 3/4. 
Next, we pick another threshold T, and we compute for each vertex v of VD(i2 U K) the excess 
= by sin. The excess measures how far the vertex deviates from the desired conflict size n/s. We 
check if we start over with a new sample. Below, we prove that 

for T = 0(s), the success probability is at least 3/4. The total success probability is 1/2, and the 
expected number of attempts is 2. Thus, in expected time 0{n log s + s log s), we can hnd a sample 

RCS with E«eVD(i?uiC) bv = 0{n) and E^,eVD(i^ui^) log 4 = 0{s). 

We now analyze the success probabilities, using the classic Clarkson-Shor method |12) . We begin 
with a variant of the Chazelle-Friedman bound m- 

Lemma 3.2. Let X he a planar point set of size m, and let T C with \Y\ < 3 and X DY = 0. 
For fixed p G (0,1], let R Y X he a random subset of size pm and let Ft! Y X he a random subset 
of size p'm, for p' = pl‘1. Suppose that p'm > 4. Fix u C X U T with |u| = 3, and let be the 
Voronoi vertex defined by u. Let bu be the number of points from XUY in the interior of the circle 
with center Vu and with the points from u on the boundary. Then, 

Pr[uu G YB{R U T)] < 646“^'’“/^ ^ YB{R' U Y)]. 

Proof. Let a = Pr[uu G VD(i?UT)] and a' = Pr[uu G VD(i?'uy)]. The vertex is in VD(i?uy) 
precisely if u C i? U T and i?u n (i? U T) =0, where are the points from X[JY inside the circle 
with center Vu and with the points from u on the boundary. If Bu H Y 0, then a = a' = 0, and 
the lemma holds. Thus, assume that B^ C X. Let du = |un W|, the number of points in u from X. 
There are ways to choose a pm-subset from X that avoids all points in Bu and contains 

all points of u n X, so 
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m — bu — du\ j ( ra 
pm — du ) I \pm 


m-bu- du- j 


du-l 

pm— du — 1 

JT pm-J 
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m — j 

j=o 

II 

o 

pm— du — 1 

W.. TT / 



m-du- j 


1 - 


j=0 


m-du-j 


Similarly, we get 


du —1 


y=n 


2=0 


pm — I 
m — i 


. p'ra—dxi—l 


n 

i=o 


1 - 


m-du- j J ' 


and since p'm > 4 and i < 2, it follows that 




/\ du p'm-du-l 


n 

i=o 


1 - 


m-du- j 


Therefore, since p' = p/2, 

du pm-du-l 


a' \p' 


n 


j=p'm-du 


1 - 


bu 


m- du- j 


, \ pm/2 

m / 


We can now bound the total expected conflict size. 


Lemma 3.3. PTe have E 


E 


veVD(i?uds') 


= Oin). 


Proof. By expanding the expectation, we get 


E 


E '>» 

DeVD(i?uds') 


^ PT[vuGyD{RUK)]bu, 

uC5Uif,|u|=3 


with Vu being the Voronoi vertex of u and bu its conflict size. By Lemma 3.2 with X 
and p = s/n, this is 


< 64 e"P ^“/2 ^ VD(i?' U K)]bu, 

uC5UiC,|u|=3 


□ 


S,Y = K 
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where -R' C 5 is a sample of size s/2. We bound this as 


i=0 uCSUi<',|u|=3 
6 ue[i,^) 

^ L D ’ 75 / 


P 


- CXJ 

< - Pr[uu G VD(i?' U K)] ^ 64e-^/^{i + 1) 

^ uCSUi4',|u|=3 j=0 

= 0{s/p) = 0{n), 


since X^uqSuA' |u|= 3 Pi’bu G VD(i?' U K)] = 0(s) is the size of VD(R' U K) and e + 1) = 

0 ( 1 ). ’ □ 


By Lemma 3.3 and Markov’s inequality, we can conclude that there is an M = 0(n) such that 


P’^[Z],;eVD(RuiC) ^ 1/4- 


Lemma 3.4. E 


I2veVD(RuK)^y^^Siv —0(s). 


Proof. By Lemma [3.2| with X = S, Y = K, and p = sjn, 


E 


tv log tv 

veVT){RUK) 


Y, Pr[uu G VD(i? U K)] t, log 

uC5Ui^,|u|=3 

< Y 64 e"P ^“/2 g YD{R' U K)]tu log t^ 

uCS'Ui^,|u|=3 

oo 

sE E 64e“t(i + l)2pr[7;^ g VD(i2'U iL)] 

i=0 uC5Ui<',|u|=3 
6ue[i,i±i) 

“ L p ’ p / 

OO 

< ^64e"t(i + 1)2 Pr[uu G VD(i2'U iL)] 

i=0 uC5UiC,|u|=3 

= 0(s). 


□ 


By Markov’s inequality and Lemma 3.4 we can conclude that there is a T = 0(s) with 


P’^[Z]i;eVD(RuiC) > 7"] < 1/4. This hnishes the first sampling phase. 

The next goal is to sample for each vertex v with > 2 a random subset Ry C By of size 
minjat^ logt-u, by} for large enough a > 0 (recall that By is the conflict set of v and that by = |i?i;|). 

Lemma 3.5. In total time O(nlogs), we can sample for each vertex v G VD(iiU K) with ty > 2 a 
random subset Ry C By of size minjat^ logt„, by}. 

Proof. First, we sample for each vertex v with ty > 2 a sequence ly of at^logt-y distinct numbers 


from {1,..., by}. For this, we use the hrst phase of the algorithm from the proof of Lemma 3.1 for 
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each such vertex, but without reusing the space. As explained in the proof of Lemma 3.1 this takes 
total time 


O 


{tv log tv ) log log ty 


o 



(4 log tv) log 


O(slogs), 


since Y2v log^?; = 0{s), and in particular tv logt^ = 0{s) for each vertex v (note that the constant 
in the 0-notation is independent of v). Also, since ty log ty = 0(s), the total space requirement 
is 0(s). 

After that, we scan through S. For each vertex v, we have a counter Cy, initialized to 0. For 
each p G S', we hnd the conflict vertices of p, and for each conflict vertex v, we increment Cy. If Cy 
appears in the corresponding set ly, we add p to Ry. The total running time is 0(nlogs), as we do 
one point location for each input point and the total conflict size is 0{n). □ 


We next show that for a fixed vertex v G VD(i? U K), with constant probability, all vertices in 
VD(i?t,) have conflict size n/s with respect to By. 


Lemma 3.6. Let v G VD(i? U K) with ty > 2, and let Ry C By be the sample for v. The expected 
number of vertices v' in VD(i?i,) with at least n/s points from By in their conflict circle is at most 
1/4. 

Proof. If Ry = By, the lemma holds, so assume that atylogty < by. Recall that ty = bys/n. We 
have 


E 


E ' 


lv'eVD{Rv) J 

b', >71/3 

v' — ' 


Y, PrKGVD(R,)], 

uCBv,\u\=3 

b[^>n/s 


where denotes the number of points from By inside the circle with center and with the points 
from u on the boundary. Using Lemma 3.2 with X = By, T = 0, and p = {atylogty)/by = 
a{s/n) logty, this is 


< Y ^ VD(R/)] 

uCBu,|u|=3 

b'y>n/s 

E PrK G VD«)] 

uCB„,|u|=3 

= 0{tf°‘/‘^ty\Ogty) < 1/4, 


for a large enough (remember that ty >2). 


□ 


By Lemma 3.6 and Markov’s inequality, the probability that all vertices from VD(R„) have at 
most n/s points from By in their conflict circles is at least 3/4. If so, we call v good, see Figure]^ 
Scanning through S, we can identify the good vertices in time O(nlogs) and space 0(s). Let s' be 
the size of VD(RU K). If we have less than s'/2 good vertices, we repeat the process. Since the 
expected number of good vertices is 35^4, the probability that there are at least s'/2 good vertices 
is at least 1/2, by Markov’s inequality. Thus, in expectation, we need to perform the sampling 
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Figure 4: A Voronoi Diagram of the sampled set R (left). The two red square vertices of VD(i?U AT) 
are not good and we need to resample within their conflict list (the blue crosses) and compute the 
new Voronoi Diagram (right). 


twice. For the remaining vertices, we repeat the process, but now we take two samples per vertex, 
decreasing the failure probability to 1/4. We repeat the process, taking in each round the maximum 
number of samples that fit into the work space. In general, if we have s'/ai active vertices in round 
i, we can take a, samples per vertex, resulting in a failure probability of 2““b Thus, the expected 
number of active vertices in round i + 1 is s'/ai^i = s'/{ai2°‘^). After 0(log* s) rounds, all vertices 
are good. To summarize: 


Lemma 3.7. In total expected time 0(n log slog* s) and space 0(s), we can find sets R ^ S and 
Rv C By for each vertex v G VD(iiU K) such that (i) |i?| = s: (ii) X]t,gvD(iJui^) ~ 

(in) for every Ry, all vertices ofYD{Ry) have at most n/s points from By in their conflict circle. 


We set i ?2 = i? U Ui;eVD(Ruii:) Lemma 3.7 |i? 2 | = 0{s). We compute VD(i ?2 U K) 


and triangulate its bounded cells. For a triangle A ot the triangulation, let r G A 2 U AT be the site 
whose cell contains A, and vi,V 2 , V 3 the vertices of A. We set B^ = {r} U lj?=i Using the next 
lemma, we show that I^IaI = 0 {n/s). 


Lemma 3.8. Let S' C and A = {ui,U 2 ,U 3 } a triangle in the triangulation o/VD(S). Let x G A. 
Then any circle C with center x that contains no points from S is covered by the conflict circles of 
vi,V 2 and U 3 . 


Proof. Let p € C and let r G S be the site whose cell contains A. We show that p is contained in 
the conflict circle of vi, V 2 , or U 3 . Consider the bisector B of p and r. Since C contains p but not 
r, we have d{x,p) < d{x,r), so x lies on the same side of B as p. Since x G A, at least one of xi, 
V 2 , X 3 , is on the same side of B as p; say vi. This means that d{vi,p) < d{vi,r), so p lies inside the 
circle around vi with r on the boundary. This is precisely the conflict circle of Xi. □ 


Lemma 3.9. Any triangle A in the triangulation 0 /VD(A 2 U K) has I^IaI = 0{n/s). 

Proof. Let x be a vertex of A. We show that by = 0{n/s). L et A r = {xi,X 2 ,X 3 } be the triangle 


in the triangulation of VD(A) that contains x. By Lemma 3.8 we have By C We 

consider the intersections By n By., for i = 1,2,3. If ty^ < 2, then = 0{n/s) and \By n Byf = 
0{n/s). Otherwise, we have sampled a set Ry^ for Xj. Let Aj = {rxi, 1 x 2 , 1 x 3 } be the triangle in the 
triangulation of VD(At,J that contains x. Again, by Lemma 3.8 we have By C By,, and thus 

is at most n/s for 

□ 


also By n By^ C ljj=i n By-. However, by construction of Ry., \By,. 0 By 
j = 1, 2, 3. Hence, \By n Byf = 0{n/s) and by = 0{n/s). 
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The following lemma enables us to compute the Voronoi diagram of R 2 D K locally for each 
triangle A in the triangulation of VD(i22 U i^) by only considering sites in B^. It is a direct 


consequence of Lemma 3.8 


Lemma 3.10. For every triangle A in the triangulation o/VD(i? 2 UiL), we have VD(S'UiL)nA = 
vd(5a) n A. 


Theorem 3 . 11 . Let S be a planar n-point set. In expeeted time 0{{n^/s) logs + n\ogs\o^ s) and 
space 0{s), we can compute all Voronoi vertices of S. 


Proof. We compute a set R 2 as above. This takes 0(n logs log* s) time and space 0{s). We 
triangulate the bounded cells of VD(i ?2 U K) and compute a point location structure for the result. 
Since there are 0{s) triangles, we can store the resulting triangulation in the workspace. Now, the 
goal is to compute simultaneously for all triangles A the Voronoi diagram VD(i?A) and to output 
all Voronoi vertices that lie in A and are defined by points from S. By Lemma |3.10[ this gives all 
Voronoi vertices of VD(5). 

Given a planar m-point set X, the algorithm by Asano et al. finds all vertices of VD(A) in 0{m) 
scans over the input, with constant workspace jl]. We can perform a simultaneous scan for all sets 
Ba by determining for each point in S all sets that contain it. This takes total time O(nlogs), 
since we need one point location for each p € S and since the total size of the Sa’s is 0{n). We 
need 0(maxA |.Ba|) = 0 {n/s) such scans, so the second part of the algorithm needs 0((n^/s) logs) 
time. □ 


As mentioned in the introduction. Theorem 3.11 also lets us report all edges of the Delaunay 
triangulation of S in the same time bound: by duality, the three sites that define a vertex of VD(S') 
also define a triangle for the Delaunay triangulation. Thus, whenever we discover a vertex of VD(S'), 
we can instead output the corresponding Delaunay edges, while using a consistent tie-breaking rule 
to make sure that every edge is reported only once. 
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A The Asano-Kirkpatrick Algorithm 


We give more details on the algorithm of Asano and Kirkpatrick jSj. Let K be a mountain with 
vertices qi,... ,qn sorted in x-order and base qiqn- We dehne the height h{qi) of g*, i = 1,...,n, as 
the distance from qi to the line through the base. Let A= (gi,..., qn) be the input array. A vertex 
qr is the nearest-smaller-right-neighbor (NSR) of a vertex qi if (i) I < r; (ii) h{qi) > h{qr)‘, and (iii) 
h{qi) < h{qk) for I < k < r. We call {qi,qr) a NSR-pair, with left endpoint qi and right endpoint 
qr- Nearest-smaller-left-neighbors (NSL) and NSL-pairs are dehned similarly. Let R be the set of 
all NSR-pairs and L be the set of all NSL pairs. Asano and Kirkpatrick show that the edges RU L 
triangulate F. We describe the algorithm for computing R. The algorithm for L is the same, but 
it reads the input in reverse. 

Let s denote the space parameter. The algorithm runs in log^ n rounds. In round i, i = 
0,... jlog^n — 1, we partition A into s* consecutive blocks of size n/sL Each block B is further 
partitioned into s consecutive sub-blocks bi,... ,bs of size In each round, we compute only 

NSR-pairs with endpoints in different sub-blocks of the same block. We handle each block B 
individually as follows. The sub-blocks of B are visited from left to right. When we visit a sub¬ 
block bj, we compute all NSR-pairs with a right endpoint in bj and a left endpoint in the sub-blocks 
6i,..., bj-i. Initially, we visit the hrst sub-block bi and we push a pointer to the rightmost element 
in bi onto a stack S. We call a sub-block with a pointer in S active. Assume now that we have 
already visited sub-blocks 6 i,..., bj-i. Let I be the topmost pointer in S, referring to an element qi 
in bji, j' < j. Furthermore, let r be a pointer to the leftmost element q^ in bj. If h{qi) > h{qr), we 
output {qi,qr) and we decrement I until we hnd the hrst element whose height is smaller than the 
current h[qi). If I leaves bj/, this sub-block becomes inactive and we remove I from S. We continue 
with the new topmost pointer as our new 1. On the other hand, if h{qi) < h{qr), we increment r 
by one. We continue until either r leaves bj or S becomes empty. Then we push a pointer to the 
rightmost element in bj onto S and proceed to the next sub-block. 

In each round, the algorithm reads the complete input once in x-order. In addition, the algorithm 
reads at most once each active sub-blocks in reverse order. Note that a sub-block becomes active 
only once. 
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